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Monte Carlo event simulation with BFKL evolution is discussed. We report current status of a 
Monte Carlo event generator ULYSSES with BFKL evolution implemented. The ULYSSES, based 
on Pythia Monte Carlo generator, would help to reveal BFKL effects at LHC energies. In particular, 
such an observable as dijet K-factor can serve as a source of BFKL dynamics at the LHC, and it 
would also help to search for new physics. 



With advent of the Large Hadron Collider (LHC) at CERN we have an opportunity to probe the Standard Model 
far beyond the explored domains. QCD is an essential ingredient of the Standard Model, and it is well tested in hard 
processes when transferred momentum is of the order of the total collision energy (Bjorken limit: Q 2 ~ s — > oc). The 
cornerstones of perturbative QCD at this kinematic regime (QCD-improved parton model)are factorization of inclusive 
hard processes [l[ and Gribov-Lipatov-Altarelli-Parisi-Dokshitzer (GLAPD) evolution equation 0]. They provide 
' a basis for the successful QCD-improved parton model. The factorization theorem [l[ for inclusive hard processes 
. , ensures that the inclusive cross section factorizes into partonic subprocess(es) and parton distribution function(s). The 
I ■ GLAPD-evolution equation governs the log Q 2 -dependence (at Q 2 — > oo) of the inclusive hard process cross-sections 
| at fixed scaling variable x ~ Q 2 / s. 

Another kinematic domain that is very important at high-energy is given by the Balitsky-Fadin-Kuraev-Lipatov 
^ \ (BFKL) limit 0,11 IE II, or QCD-Regge limit, whereby at fixed Q 2 > A^ CD , s -> oo. In the BFKL limit, the BFKL 
II ■ evolution in the leading-log approximation (LLA) governs log(l/a;) evolution (at x — > 0) of inclusive processes. 
Qh' One of the key BFKL features is a relation of the highest eigenvalue, w mm , of the BFKL equation (EEL IE [with 
= to] the intercept of the Pomeron, which in turn governs the high-energy asymptotics of the total cross-sections: a ~ 
r"" 1 1 ' (s/so) Qjp_1 = (s/s ) w " ax , where the Regge parameter so defines the approach to the asymptotic regime. The BFKL 
Pomeron intercept in the LLA turns out to be rather large: aip — 1 = u) T? x a = 12 log 2 (as/ir) — 0.54 for as = 0.2; 
hence, it is very important to take into account NLLA corrections M, |^, I^Jld. ITlT| to the BFKL. Note that the BFKL 
evolution in the next-to- leading-log approximation (NLLA) 0, H, unlike the LLA BFKL 0,0,0, includes GLAPD 
evolution with the running coupling constant of the leading-order (LO) GLAPD, as(Q 2 ) = 47r//3 log(Q 2 /Ag CI? ). 
One of the striking features of the NLLA BFKL analysis 0] is that the NLLA value for the intercept of the 
H ' BFKL Pomeron, improved by the BLM procedure [HI, has a very weak dependence on the gluon virtuality Q 2 : 
ajp — 1 = wJvlla — 0-13 - 0.18 at Q 2 = 1 - 100 GeV 2 . This agrees with the conventional Regge theory where 
one expects universal intercept of the Pomeron without any <5 2 -dependence. The value of NNLA BFKL Pomeron 
intercept becomes compatible with the available data on hard QCD Pomeron. So, NLLA BFKL approach 0, [l2j] 
posses all basic features of BFKL, but includes also running coupling effects and realistic value of the hard Pomeron 
IJ • intercept. 

. i-H ! Therefore, the BFKL and especially the NLLA BFKL 0, H, [12] are anticipated to be important tools for exploring 
■ the high-energy limit of QCD. 

It should be stressed that in contrast to the GLAPD, BFKL dynamics involves parton distributions unintegrated 
over k t - So, to reveal BFKL effects one needs to deal, e.g., with parton (jet) production. For jet production in the 
LLA BFKL within fc t -factorization [3] one can use effective Feynman-like rules for inclusive jet cross sections [is! ]. 

However, extremely sophisticated design of contemporary high-energy detectors requires from theory and phe- 
nomenology a Monte Carlo event generator implementation. Widely used Monte Carlo event generator Pythia (l6j 
contains basically GLAPD-evolution. To incorporate BFKL-evolution Pythia has been modified in parton shower 
kernel and parton distribution functions. 

For Monte Carlo event simulations of GLAPD evolution it is convenient to use the Dokshitzer-Diakonov-Troyan 
(DDT) [13 representation: 
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where P a b{z) arc kernels (splitting functions) of GLAPD equation for partons a and b, 

rn 2 



a(x,(J, ) 



dk t f a (x,k t ,fx ), a = g,q,q 



is GLAPD parton distribution for parton a, f a {x,k 2 ,^i 2 ) is so-called unintegrated parton distribution used in BFKL 
approach and 



T a {klf) ^ exp (- jf dn 2 ^1 Yl J/t CA-(0 



is QCD analog for Sudakov form factor (DDT form factor which defines the probability not to emit a gluon. 

In terms of unintegrated parton distribution the GLAPD equation can be represented as: 
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Then, a unified BFKL-GLAPD equation can be presented in the following way [18| : 



/ ff (z,fc 2 ,/i 2 ) = T g (k u n) 
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The Eq. (fTJ) is similar to Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [T^| equation, which is a version of BFKL 
equation with imposed angle ordering. 



Inclusive to exclusive ratio 
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FIG. 1: Dijet K-factor at LHC: inclusive dijet cross section over only-two-jet cross section ratio. The results are shown with 
full CMS detector simulation (arbitrary normalization) for luminosity 20 pb" 1 . 
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Monte Carlo event generator Ulysses [2(| contains an implementation of the Eq. |l| into Pythia [l6[ parton shower 
simulation and embedded BFKL unintegrated parton distributions. Preliminary results for dijet K-factor at CMS 
detector of the LHC are shown in Fig. 1, where dijet K-factor for jets with E t > E™ m = 60 GeV is defined as ratio 
of inclusive dijet cross section to "exclusive" one. "Exclusive" dijet event means the event when only two jets with 
E t above E™ m are available ("Born" cross section). 

To summarize, BFKL Monte Carlo should help to look at LHC for BFKL dynamics as a new feature of asymptotic 
QCD. It also has a potential for search of new physics, e.g., graviton production in trans-Planckian regime [2l|. 
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